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ASTEROSEISMIC CONSTRAINTS FOR GAIA 
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Abstract. Distances from the Gaia mission will no doubt improve our understanding of stellar physics 
by providing an excellent constraint on the luminosity of the star. However, it is also clear that high 
precision stellar properties from, for example, asteroseismology, will also provide a needed input constraint 
in order to calibrate the methods that Gaia will use, e.g. stellar models or GSP_Phot. For solar-like stars 
(F, G, K IV/V), asteroseismic data delivers at the least two very important quantities: (1) the average 
large frequency separation (Au) and (2) the frequency corresponding to the maximum of the modulated- 
amplitude spectrum i/ mlx . Both of these quantities are related directly to stellar parameters (radius and 
mass) and in particular their combination (gravity and density). We show how the precision in (Au), f max , 
and atmospheric parameters T e s and [Fe/H] affect the determination of gravity (log<?) for a sample of well- 
known stars. We find that \ogg can be determined within less than 0.02 dex accuracy for our sample while 
considering precisions in the data expected for V ~ 12 stars from Kepler data. We also derive masses and 
radii which are accurate to within la of the accepted values. This study validates the subsequent use of 
all of the available asteroseismic data on solar-like stars from the Kepler field (> 500 IV/V stars) in order 
to provide a very important constraint for Gaia calibration of GSP_Phot through the use of logg. We note 
that while we concentrate on IV/V stars, both the CoRoT and Kepler fields contain asteroseismic data on 
thousands of giant stars which will also provide useful calibration measures. 
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1 Introduction 

The ESA Gaiffl mission is due to launch in Autumn 2013. Its primary objective is to perform a 6-D mapping of 
the Galaxy by observing over 1 billion stars down to a magnitude of V = 20. The mission will yield distances to 
these stars, and for about 20/100 million stars, distances with precisions of less than 1%/10% will be obtained. 

Gaia will obtain its astrometry by using broad band "G" photometry (similar to a V magnitude). The 
spacecraft is also equipped with a spectrophotometer comprising both a blue and a red prism BP/RP, delivering 
colour information. A spectrometer will be used to determine the radial velocities of objects as far as G = 17 
(precisions from 1-20 kms" 1 ), and for stars with G < 11 high resolution spectra (R~ll,500) will be available. 

One of the main workpackages devoted to source characterisation is GSP_Phot whose objectives are to obtain 
stellar properties for 1 billion singl e stars by using t he G band photometry, the parallax n, and the spec- 
trophotometric information BP /RP ()Bailer-Jonesll201oT) . The stellar properties that will be de rived are e f fective 
temperature T c g, extinction Aq in the G band, surface gravity logg, and metallicity [Fe/H] . ILiu et al. (2012) 



compare different methods to determine these parameters and they estimate typical precisions in log g on the 
order of 0.1 - 0.2 dex for main sequence late-type stars, and mean absolute residuals (true value minus inferred 
value from simulations) no less than 0.1 dex for stars of all magnitudes. 

A calibration plan using forty bright benchmark stars has been put in place to deliver the best stellar 
models. These will be used on ~5000 calibrations stars which will be observed by Gaia. However, for most of 
these fainter stars logg remains quite unconstrained, and this will inherently reduce the full capacity of source 
characterisation with Gaia data. 

In the last decade or so, much progress in the field of observational asteroseismology has been made, especially 
for stars exhibiting Sun-like oscillations. These stars have deep outer convective envelopes where stochastic 
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turbulence gives rise to a broad spectrum of excited resonant oscillation modes e.g. Brown fc Gillilaiidl 19941 
The power spectra of such stars can be characterised by two mean seismic quantities: (Au) and i/ max . The 
quantity (Av) is the where z^„ is a frequency 

with degree / and radial order n, and v maK is the frequency corresponding to the maximum amplitude of 
the bell-shaped frequency spectrum. The following scaling relations have also been shown to hold: Eq. 1: 
(Au) w M°- 5 J?- 15 (Az/) s and Eq. 2: iw » MRr *{T eS /5777)- a5 tw.^ (Eq. 1) where (Ai/) Q = 134.9 /jRz 
and y maX:Q = 3,050 /xHz ( Kieldsen fc Bedding|ll995h . 

Of particular interest for Gaia is the Kepler (http : //kepler .nasa. gov) field of view — ~100 square-degrees, 
centered on galactic coordinates 76.32°, +13.5°. Kepler is a NASA mission dedicated to characterising planet- 
habitability ( Borucki et al. 2010l ). It obtains photometric data of ~150, 000 stars with a typical cadence of 30 
minutes. A subset of stars (< 1000 every month) acquire data with a point e very 1 min ute. This is sufficient 
to detect and characterise Sun-like oscillations in many stars. Verner et al. ( 2011 ) and IChaplin et aL ( 2011 ) 
recently showed the detections of these mean seismic quantities for a sample of >500 F, G, K IV/V stars with 
typical magnitudes 7 < V < 12, while both CoRoT and Kepler have both shown their capabilities of detectin g 
these same seismic quantities in 1000s of red giants (jHekker et al. 1 120091: iBaudin et al1l201lt iMosser et al]|2012h . 

With the detection of mean seismic quantities in hundreds of stars, the Kepler field is very promising for 
helping to calibrate the GSP_Phot methods. In particular, t hey deliver one o f the four stellar properties to be 
extracted by automatic analyses of Gaia data, namely log g. iGai et al.1 (|2011l ) studied the distribution of errors 
for a sample of simulated stars using seismic data and a grid-based method based on stellar evolution models. 
They concluded that a seismic log g is almost fully i ndepe ndent of the input physics in the stellar evolution 
models that are used. More rece ntlv lMorel fc Migliol (|2012f) compared classical determinations of log g to those 
derived alone from the scaling relation (Eq. [2]), and concluded that the mean differences between the various 
methods used is ^0.05 dex, thus supporting the validity of a seismic determination of log g. However, to date, 
no study has been done to validate the accuracy of a seismic log g (how closely it resembles the true value) by 
using stars with measured radii and masses. This is the objective of this work. 



2 A comparison of the direct and seismic methods for determining \ogg. 
2.1 Observations and direct determination of log g 

We aim to compare an asteroseismically derived log<? with the true known value for a sample of stars. We chose 
a sample of seven bright well-characterised stars for which the radius is known via interferometry or a binary 
solution and the mass is known from either the binary solution or a detailed seismic analysis. Table [T] lists the 
sample of stars along with the observed values of (Au), u max , T e ff, [Fe/H], M, and R. The final column in the 
table gives the true value of logg derived from M and R. 



2.2 Seismic method to determine log g 



We use a grid-based method, RadExlO, to determine an astero seismic value of logq (ICreevev et al.ll2012T) . The 
grid was constructed using the ASTEC stellar ev olution code ( Christensen-Dalsgaard 20081 ) without diffusion 
effects and the same input physics as described in Creevev et al.l ( 2012 ). 

The grid considers models with masses M from 0.75 - 2.0 M in steps of 0.05 M , ages t from ZAMS to 
subgiant, the initial chemical composition Z\ (metallicity) spans 0.007 - 0.027 in steps of ~ 0.003, while X\ 
(hydrogen) is set to 0.70: this corresponds to an initial He abundance Y\ = 0.263 — 0.283. The mixing length 
parameter a — 2.0 is used, which was obtained by calibrating it with solar data. 

To obtain the grid-based model stellar properties (log g, M, R, L, t) we perturb the set of input observations 
using a random Gaussian distribution, and compare the perturbed observations to the model ones. The input 
observations consist primarily of (Au), i/ max , T Q g, and [Fe/H], although other inputs are possible, for example, 
L or R. The stellar parameters and uncertainties are defined as the mean value of the fitted parameter from 
10,000 realizations, with the standard deviations defining the la uncertainties. 



2.3 Analysis approach 

We determine a seismic log<? for the stars using the method explained above, and the following data sets: 

(51) {(Ai/), iw, T off , [Fe/H]}, 

(52) {(Au), z/ max , T G ff}, and 
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Table 1. Observed properties of the reference stars 



Star 


(Av) 


^max 


T c s 


[Fe/H] 


R 


M 


log.g 




(jjlHz) 


(mHz) 


(K) 


(dex) 


(R®) 


(M ) 


(dex) 


aCenB 


161.5±0.11 la 


4.0 lQ 


5316±28 16 


0.25±0.04 16 


0.863±0.005 lc 


0.934±0.0061 ld 


4.538±0.008 


18 Sco 


134.4±0.3 2a 


3.1 2q 


5813±21 2a 


0.04±0.01 2a 


1.010±0.009 2a 


1.02±0.03 2a 


4.438±0.005 


Sun 


134.9±0.1 3a 


3.05 3b 


5778±20 3c 


0.00±0.01 3d 


1.000±1.010 3d 


1.000±0.010 3d 


4.438±0.002 


aCenA 


105. 6 4a 


2.3 4q 


5847±27 lfc 


0.24±0.03 lb 


1.224±0.003 lc 


1.105±0.007 ld 


4.307±0.005 


HD 49933 


85.66±0.18 5a 


1.8 5q 


6500±75 56 


-0.35±0.10 5fc 


1.42±0.04 5c 


1.20±0.08 5c 


4.212±0.039 


Procyon 


55.5±0.5 6a 


1.0 6b 


6530±90 6c 


-0.05±0.03 6d 


2.067±0.028 6e 


1.497±0.037 6/ 


3.982±0.016 


/3Hydri 


57.24±0.16 7a 


1.0 7a 


5872±44 76 


-0.10±0.07 7c 


1.814±0.017 7h 


1.07±0.03 7b 


3.950±0.015 



References: HKieldsen et al l l|2005l ). ^Porto de Mello et al l d200Sh. 1 1Kervella et all (|2003h. ^ourbaix et all (|2002l ), 
"iBazot et al.1 (|201ll) , 3a Ta king the average of Table 3 from iToutain fc Froehlichl l| 19921), 3 lKieldsen & Beddind (|l995h . 
3 lGrevesse fc Sauvall (|l998l ), 3d We adopt a typical error o f 0.01 in [Fe/H] M and R 4 °lBouchv fc Carrier! (|2002T l. 5a Using 
the I = modes wit h Height / Noise >1 fro m Table 1 of iBenomar et al.l l|2009l), 5 lKallinger et all l|2010l) Z = 008 ± 
0. 002 is referenced, 5t u3igot et al.1 l|201lf). 6t iEggenberger et al.1 (|2004l), 6 lMartic et all (|2004t), 6 lFuhrmann et alJ dl997h. 
lAllende Prieto et al l <l2002h . 6 lKervella et~all (120041 1. 6J lGirard et all ^OOd l. 7 ° lBedding et~ai1 (|2007h . 7 lNorth et all 
20071 ). 7 lBruntt et al] (|2010h . 
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(S3) {(Ai/), ^ max }. 

For the potential sample of Gaia calibration stars, [Fe/H] is not always available, and in some cases, a photo- 
metric T e ff may have various estimations. For these reasons we include S2 and S3. 

The observational errors in our sample are very small due to the brightness and proximity of the star, so 
we also derive an a steroseismic log g while considering observational errors that we expect for Kepler stars (see 
IVerner et ai1l201lh . We consider three types of observational errors: 
(El) the true measurement errors from the literature, 

(E2) typically "good" errors, i.e. cr((A^),j/ max ,r cff ,[Fe/H]) = 0.5 /jHz, 5%, 70 K, 0.08 dex, 

(E3) "not-so-good" errors (e.g. V-11,12), a((A^)Vmax,^eff,[Fe/H]) = 1.1 /iHz, 8%, 110 K, 0.12 dex. 



2.4 Seismic versus direct log<? 

In Figure [T] we compare the asteroseismic logg with the true logg for the seven stars. Each star is represented 
by a point on the abscissa, and the y-axis shows (seismic - true) value of logg. There are three panels which 
represent the results using the three different subsets of input data. We also show for each star in each panel 
three results; in the bottom left corner these are marked by 'El', 'E2', and 'E3', and represent the results using 
the different errors in the observations. The black dotted lines represent (seismic - true) logg = 0, and the grey 
dotted lines indicate ±0.01 dex. 

Figure [1] shows that for all observational sets and errors log g is generally estimated to within 0.02 dex in both 
precision and accuracy. This result clearly shows the validity of the mean seismic quantities and atmospheric 
parameters for providing an extremely precise value of logg. Other general trends that can be seen are that 
the typical precision in logg decreases as (1) the observational errors increase (from El - E3), and (2) the 
information content decreases (SI - S2 - S3, for example). One noticeable result is the systematic offset in the 
derivation of log<? for HD 49933 when we use [Fe/H] as input (SI). This could be due to an incorrect metallicity, 
an error in the adopted true logg or a shortcoming of the grid of models. 

Figure [2] shows the seismic radius and mass determinations of the sample stars using SI while considering the 
three sets of errors. We find that with good observational errors, the radii are matched to within 1% (accuracy) 
with typical precisions of 2-3%, while the masses are matched to within 1-4% with typical precisions of 4-7%. 
Here it can be seen that the offset found for HD 49933 in log g is related to the reference mass value (the 
radius seems to be consistent). For S2 and S3 the uncertainties begin to grow very large; 2-5% and 3-10%, 
respectively in radius, and 5-15% and 10-35% in mass, while the accuracies also decrease (not shown) although 
to within < 1.5c for all results. These results indicate that for the most precise determination of mass and 
radius, a seismic index and both T c g and [Fe/H] are necessary, unlike log g where the seismic information alone 
(or including T e ff) can produce an accurate result. 
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Fig. 1. Seismic-minus-true log g for the seven sample stars while considering different sets of input observations (different 
panels) and different observational errors (El, E2, E3). 



Table 2. Stellar properties for the reference stars derived by RadExlO 

Star logg(dex) R (R Q ) M (M Q ) L (L ) Age (Gyr) 

a Cen B 4.527 ± 0.004 0.859 ± 0.007 0.905 ± 0.023 0.52 ± 0.02 9.4 ± 2.0 

18 Sco 4.441 ± 0.004 1.018 ± 0.008 1.042 ± 0.019 1.07 ± 0.04 4.8 ± 0.9 

Sun 4.438 ± 0.001 1.000 ± 0.002 1.000 ± 0.005 1.01 ± 0.03 6.3 ± 0.6 

a Cen A 4.312 ± 0.004 1.223 ± 0.010 1.119 ± 0.024 1.56 ± 0.08 7.0 ± 0.9 

HD 49933 4.195 ± 0.007 1.418 ± 0.022 1.148 ± 0.054 3.23 ± 0.22 3.5 ± 0.6 

Procyon 3.981 ± 0.006 2.072 ± 0.024 1.497 ± 0.041 7.08 ± 0.54 2.1 ± 0.2 

/3 Hydri 3.957 ± 0.010 1.840 ± 0.045 1.119 ± 0.086 3.55 ± 0.31 6.8 ± 1.0 



2.5 Systematic errors in observations 



To study the effect of systematic errors in the atmospheric parameters, we repeated our analysis for /3 Hydri 
using th ree sets of in p ut da ta that change only in T e g and [Fe/H] considering the E2 errors. The first set (1) 
uses the iNorth et all (120071) values (5872, -0.10), the second set (2) uses (5964, -0.10), and the third set (3) 



uses 



da Silva et al 



(|2006l) values (5964, -0.03). For SI we derived logg = 3.96, 3.97, and 3.97 dex for case 1, 
2, and 3, respectively (ref. value is 3.95 dex). Excluding [Fe/H] (S2) we derived logg = 3.95 and 3.96 for case 
1 and 2, respectively. Here we can conclude that errors in the atmospheric parameters can change log g by up 
to 0.02 dex, and in the absence of an accurate [Fc/H] it is better to exclude it. To determine the mass, radius, 
and age, however, [Fe/H] is a very imporant constraint. 



4 - 
2 - 




aCenB 18 Sco Sun aCenA HD49933Procyon /SHydri 



aCenB 18 Sco Sun aCenA HD49933Procyon (3Hydri 



Fig. 2. Seismic-minus- True values of radius (left panel) and mass (right panel) using SI. 
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3 Conclusions 

We summarize the stellar properties of the sample stars in Table [2] derived by RadExlO using (Ais), ^ max , 
T e g, and [Fe/H], and the true observational errors. We highlight the excellent agreement between seismically 
determined parameters and those obtained by direct mass and radius estimates (compare Tables Q] and [2]) . In 
only two cases (a Cen B and 18 Sco), we find that log<? and mass are determined with a difference of just 
over la for SI and S3, while for S2 we find that logg is accurate to within its a for all stars. This study 
validates the accuracy of a seismically determined log g while also highlighting the excellent precision that can 
be obtained using seismic data. If we relax the observational errors to those typical of what is available for the 
sample of ^500 Kepler F, G, K IV/V Kepler stars, then we obtain \ogg with precisions of less than 0.02 dex 
for SI (including [Fe/H] as a measurement) and less than 0.03 dex for S2 (excluding [Fe/H]) for even "poor" 
observational errors on the input seismic and atmospheric data. We also showed that we can expect to find a 
typical systematic error of no bigger than 0.02 dex arising from an error in the atmospheric parameters. 

OLC is a Henri Poincare Fellow at OCA, funded by the Conseil General des Alpcs-Maritimes and OCA. 
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